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Nonlinear Growth of Firehose and Mirror Fluctuations in Astrophysical Plasmas 
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In turbulent high-beta astrophysical plasmas (exemplified by the galaxy cluster plasmas), 
pressure-anisotropy-driven firehose and mirror fluctuations grow nonlinearly to large amplitudes, 
8B/B ~ 1, on a timescale comparable to the turnover time of the turbulent motions. The principle 
of their nonlinear evolution is to generate secularly growing small-scale magnetic fluctuations that 
on average cancel the temporal change in the large-scale magnetic field responsible for the pres- 
sure anisotropies. The presence of small-scale magnetic fluctuations may dramatically affect the 
transport properties and, thereby, the large-scale dynamics of the high-beta astrophysical plasmas. 

PACS numbers: 52.35.Py, 94.05.-a, 95.30.Qd, 98.65.Hb 



Introduction. — Many astrophysical plasmas are mag- 
netized and weakly collisional, i.e., the cyclotron fre- 
quency flj is much larger than the collision frequency 
Vn and the Larmor radius pi is smaller than the mean 
free path A m f p . In such plasmas, all transport properties, 
most importantly the viscosity and thermal conductivity, 
become anisotropic with respect to the local direction of 
the magnetic field [l[ — even if the field is dynamically 
weak. 

As a typical example where just such a physical situ- 
ation is present is galaxy clusters @, 0] . While param- 
eters vary significantly both within each cluster and be- 
tween clusters, the weakly collisional magnetized nature 
of the intracluster medium (ICM) is well illustrated by 
the core of the Hydra A cluster, where ~ 1CP 2 s _1 , 
u ti ~ 10" 12 s" 1 and pi ~ 10 5 km, A mfp - 10 15 km 0]. 
Modeling global properties of clusters and physical pro- 
cesses inside them, such as shocks, fronts, radiobubbles, 
or the heating of the ICM [f|, can only be successful if 
the viscosity and thermal conductivity of the ICM are 
understood [f| . Another fundamental problem is the ori- 
gin, spatial structure and the global dynamical role of 
the magnetic fields in clusters. Turbulent dynamo mod- 
els again require knowledge of the ICM viscosity [!, 0, 0] , 
which itself depends on the field structure, so the prob- 
lem is highly nonlinear and is as yet unsolved. 

An additional complication is that in a turbulent 
plasma, pressure anisotropies develop in a spontaneous 
way [2, [1, HI- In high-beta plasmas, they trigger a num- 
ber of instabilities, most interestingly, firehose and mir- 
ror 0, [l3| ■ The instabilities are very fast compared to 
the motions of the ICM and give rise to magnetic fluc- 
tuations at scales as small as pi. The spatial structure 
and the saturated amplitude of these fluctuations must 
be understood before quantitative models of transport 
can be constructed. In this Letter, we demonstrate how 
the nonlinear kinetic theory of these fluctuations can be 
constructed, elucidate the basic physical principle behind 



their nonlinear evolution and show that they do not satu- 
rate at small quasilinear levels [IH , but grow nonlinearly 
to large amplitudes (5B/B ~ 1). 

The physical origin of pressure anisotropies. — A fun- 
damental property of a magnetized plasma is the conser- 
vation of the first adiabatic invariant for each particle, 
p = v 2 ± /2B (on time scales ^> ft" 1 ). This implies that 
any change in the field strength must be accompanied by 
a corresponding change in the perpendicular pressure, 
Pa_/ B const. In a heuristic way, we may write Q 
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where the last term represents collisions relaxing the 
pressure anisotropy. On the other hand, the magnetic 
field is frozen into the plasma flow velocity u and the 
field strength obeys [l2j 
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where d/dt = d/dt + u ■ V, b = B/_B and 70 is the 
turnover rate of the turbulent motions. Taking the two 
terms in the right-hand side of Eq. to be comparable 
and using Eq. (2j), we get A = (pj_ -p\\)/p± ~ Jo/fa- 
This is the anisotropy persistently driven by the turbu- 
lent motions, which are excited at the large (system-size) 
scales by various macroscopic mechanisms [7] . 

If the turbulence is Kolmogorov, the dominant contri- 
bution to the turbulent stretching and, therefore, to the 
pressure anisotropy, comes from the viscous scale l v ~ 
Re~ 3 / 4 L, where L is the outer scale. The viscous-scale 
motions have the characteristic velocity u ~ Re -1 / 4 £/, 
where U is the characteristic velocity at the outer scale. 
The Reynolds number Re = UL/v is calculated using the 
viscosity of an unmagnetized plasma v ~ fthiA m f p (wthi 
is the ion thermal speed) because for the motions that 
change the field strength the viscosity is not reduced by 
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the magnetic field 0, 0] . We now introduce a small pa- 
rameter e ~ MRe -1 / 4 , where M = U/v^i is the Mach 
number 0. Then we can order u/vthi ~ A m f p /L ~ e, 
whence A — To/f-H ~ uA m f p /Z„Vthi ~ e 2 - Using again 
our fiducial parameters for the Hydra A cluster core, 
U ~ 250 km/s, v tu ~ 700 km/s, L - 10 17 km 0, 
we have e ~ 0.1. The relatively small resulting typical 
anisotropy due to turbulence will have a dramatic effect 
on the magnetic field. 

Qualitative derivation. — Consider first the firehose in- 
stability. It is activated when, or in the regions where, 
A < 0, i.e., the magnetic-field strength is decreas- 
ing. Such events/regions will always exist in a turbulent 
plasma. The growing fluctuations are polarized as Alfven 
waves, with magnetic perturbations perpendicular to the 
original field: B = B + (JB_l. Using Eq. Jl}, we estimate 
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where 70 = (l/Bo)dBo/dt < 0, the instability growth 
rate is 7 = (|A| - 2//? 4 ) 1 / 2 fc j |V thi > 7o (for k Pl < 1) 
Pi = ^ nm i nv thi/Bo> an d the overbar denotes aver- 
aging over the fluctuation scales. Intuitively, the fluctua- 
tions are averaged because the particles streaming along 
the field lines traverse the field fluctuations faster than 
the fluctuations grow (k\\Vthi 7)- Initially, 7 ^> vu\ 
as SB± grows, the instability is quenched because the 
negative anisotropy associated with the large-scale tur- 
bulence is compensated by a positive anisotropy due to 
the small-scale fluctuations. The amplitude at which the 
quenching occurs is SBj_/B ~ (Ito 1/ ^ii) 1 ^ 2 ~ e. This 
estimate can also be obtained via a formal quasilinear 
calculation 11]. However, it does not, in fact, describe 
a steady state. Indeed, if 5B± stops changing while 
the unperturbed field B continues to decrease, the re- 
sulting negative pressure anisotropy is again uncompen- 
sated and the firehose instability will be reignited. Since 
the anisotropy is reduced in the nonlinear regime, the 
growth of the fluctuations eventually slows down so that 
7 <C vu. Then Eq. ([3]) shows that the anisotropy stays 
at the marginal level if (1/ ' B^dSB 2 ^/ 'dt ~ I'yo ] 5 whence 
5B±/Bq ~ (Ito I*) 17 ' 2 ■ The physical principle of this non- 
linear evolution is that the to tal average field strength 
does not change: d(B$ + SB±)/dt = 0. 

Thus, after an initial burst of exponential growth, the 
firehose fluctuations grow secularly until the anisotropy- 
driving fluid motion decorrelates. As this happens on 
the time scale ~ |7o| -1 > the fluctuations will have time 
to become large, 5B±/Bo ~ 1. For Hydra A parameters 
used above, the time needed for that is |7o| _1 ~ 10 6 yrs. 

Kinetic theory. — We now derive these results in a sys- 
tematic way. Although finite ion Larmor radius (FLR) 
effects are important for the quantitative theory of the 
firehose instability [3 EH' the limit kpi <§; 1 provides the 
simplest possible analytical framework for elucidating the 



key elements of the nonlinear physics, which persist with 
FLR 0. We start with the Kinetic MHD equations [13], 
valid for kpi <C 1 and u> <C fi»: 
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B Vu. 
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(5) 



We set n = const and V • u = 0. This can be ob- 
tained self-consistently, but to reduce the amount of 
formal derivations we simply assume incompressibility 
at all scales (the motions are subsonic). The pressure 
anisotropy is p± — P\\ = J d 3 v m i (v 2 _/2 — v 2 )f(t,r,v), 
where / is the ion distribution function and v the ion ve- 
locity in the frame moving with the mean velocity u. The 
electron contribution to J>ii is smaller by (me/m^) 1 / 2 . 
The ion distribution function satisfies [171 ] 

dt rrii J \ dv v d£ 



bb : Vu 



l-3£ 2 df 3, ~ df 



,(6) 



where v = |v|, £ = v\\/v and the last term on the left- 
hand side is the collision operator. 

We take B = B + <$Bj_, u = u + 8u±, and £ji = 0, 
where the slow fields Bo , uo (the background turbulence) 
vary at the rate 70 on the scale l v of the viscous motions 
(or larger) and the fast perturbations <5B_l, Su± have the 
growth rate 7 and wavenumber k. We formally order 
all scales and amplitudes with respect to the small pa- 
rameter e introduced above. As we see from Eq. ©, it is 
sensible to let the fluctuation growth rate be (at least) the 
same order as the collision rate: 7 ~ ekvthi ~ vu, whence 
fcii ~ (eAmfp) -1 [18]. For the fluid motions, Uo/^thi ~ e, 
70 ~ e 3 fcw t hi, and l^ 1 ~ e 2 k. The expected fluctua- 
tion level at which the instability starts being nonlinearly 
quenched tells us to order SB±/Bq ~ e and, using Eq. 
([5]), 5u±/v t hi ~ e 2 - From Eq. (gj) we see that the pressure 
anisotropy is destabilizing only if it is not overwhelmed 
by the magnetic tension, so we order l//3j ~ A ~ e 2 . 

We seek the distribution function / = fo + Sf + 
Sf + • • • , where /q only has slow variation in space and 
time. To order e (the lowest nontrivial order), Eq. ^ 
becomes £ubo • V<J/i — (dfo/dt) c — 0. Averaging along 
the magnetic field, we get (dfo/dt) c = 0, whence fo is 
a Maxwellian: f = n exp(— v 2 /v 2 hl )/ (Trv 2 hi ) 3 / 2 . Then 
£vbo • V<5/i = 0, i.e., Sf has no fast variation along the 
magnetic field. To order e 2 , we learn, in a similar fashion, 
that Sf converges to a Maxwellian on the collision time 
scale (so it can be absorbed into fo) and that Sf has no 
fast variation along the magnetic field. Finally, to order 
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e 3 , the kinetic equation averaged along bp is 



dSf 2 fdSf 2 



dt 



V dt 



= bb : Vu (1 - 3£ 2 



V f 
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where the overbar denotes spatial averaging along the 
field line. In order to solve this equation, we assume, 
as a simple model, a pitch-angle-scattering collision op- 
erator with a constant collision rate: (dSf2/dt) c = 
(v u /2)(d/d£){l - ( 2 )d8f 2 /d£,. While this is not quanti- 
tatively correct, it is sufficient for our purposes. Solving 
for 6/2 and calculating the pressure anisotropy, we find 
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= 3 f dt'e-^^-^bb : Vu(t') 
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where 70 = (1/ Bo)dBo / dt < 0, pa = nomiV^ hi /2, and we 
have used Eq. ([2]) . Eq. {§]) is the quantitative form of Eq. 
<j3j> . Note that it generalizes the Braginskii |l| formula 
p± — pii = (p /vii)hb : Vu, which is only valid for fields 
varying slowly in space in time (cf. Q). 

Applying our ordering to Eqs. (U|) and $5§, we get 



d 2 6B A 



A(t) 
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(9) 



where /3j = 8ttpq/B 2 and V|| is the gradient along Bo. 

Eqs. ((9|) and ([8]) describe the evolution of firehose 
perturbations both in the linear and nonlinear regimes. 
Consider the evolution of a single Fourier mode (Fig.[T]). 
When SB^/Bq < e, the first (linear) term in Eq. © 
dominates and the perturbations grow exponentially with 
the firehose growth rate 7 = (|7o|/^ii — 2//? i ) 1 / 2 fc||i;thi- 
Once the nonlinearity becomes significant (at 
5B±/Bq ~ e), the anisotropy is gradually suppressed 
and, for tv u > 1, SB 2 ± {t)/B^ ~ At - 2v u / 'k 2 v 2 hi t + 



ci 



At J dt 



e -3v«i72 t /-l/4 



sm 



% /2A7 7 3A;||U thi t' 3 / 2 + c 2 ^ 

where A = 2(|7o| — 2vu/ @i) and Ci and C2 are integration 
constants. The dominant behavior (the first term) is 
the secular growth we already derived qualitatively 
above. The second term is the long-time subdominant 
correction and the third is an oscillatory transient, which 
decays on the collision time scale. 

Considering the nonlinear evolution from arbitrary 
initial conditions involving many Fourier modes requires 
inclusion of the FLR terms that set the wave number 
of maximum growth. While the spatial structure of the 
fluctuations becomes more complex and a power-law 
energy spectrum emerges (la |. the key physical result 
derived above persists: the fluctuation energy grows 
secularly with time until finite amplitudes are reached. 

The mirror instability. — The nonlinear evolution of the 
mirror instability shares some of the features of the fire- 
hose, but the full kinetic calculation is much more com- 
plicated. Here we only present a qualitative discussion. 
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FIG. 1: Evolution of 5Bj_(t)/Bo and A(t) (inset) obtained 
by numerically solving Eqs. (J5J) and (JHJ for a single Fourier 
mode. Here 'yo/va = 0.01, f3i = 1000 and k\\vthi/V2 Vu = 10. 

The mirror instability is triggered for A > (increas- 
ing B), has the growth rate 7 ~ Afeii^thi for kpi -C 1, 
and gives rise to growi ng p erturbations of the magnetic- 
field strength, 5Bu 0, [h|. The pressure anisotropy is, 
as before, determined by the changing field strength seen 
on the average by parallel-streaming particles: 



7o 



1 



d SB, 



7 + vu dt B 



(10) 



For particles traveling the full length of the field line, 
5B^ = 0; the particles for which £ < £ tr - \SB\\/ B \ 1 / 2 
are trapped by the fluctuations ("mirrors") and play a 
key role in the nonlinear dynamics [19]. Trapping be- 
comes important when the bounce frequency approaches 
the instability growth rate: uj^ ~ ^11 ^thiCtr ~ 7 ~ 
(7oM<)A||Vthij or 5B\\/B Q ~ {lo/va) 2 ~ e 4 . For 
amplitudes above this level, SBh/Bq ~ ^tiSB^/Bo ~ 
— |<5_B||/i?o| 3 / 2 (negative because particles are trapped in 
the regions of weaker field). We substitute this estimate 
into Eq. (|10[) . assume slow evolution (7 <C Vu), and 
find that the marginal state is achieved for 5B\\/Bq ~ 
(lot) 2 / 3 - This secular growth continues until the turbu- 
lent motion responsible for the pressure anisotropy decor- 
relates, by which time 8B\\/Bq ~ 1. The FLR effects, 
while important jnj [2(| , are ignored in this qualitative 
argument, but are unlikely to change the main result (sec- 
ular growth). 

Conclusion. — We have shown that, in high-beta turbu- 
lent plasmas, small-scale magnetic fluctuations are con- 
tinually generated by plasma instabilities and grow non- 
linearly to large amplitudes, SB/B ~ 1, so strongly 
"wrinkled" magnetic structures emerge on the fluid time 
scales. The main difference between our theory and most 
others 11, 15, [ijj [2(| is that they consider an initial 
pressure anisotropy gradually cancelled by fluctuations 
in a collisionlcss plasma, whereas in our calculation, the 
anisotropy is continually driven by the turbulent motions 
and relaxed by (weak) collisions; the evolution of the fluc- 
tuations is followed over times longer than the collision 
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time, up to the fluid time scale. The underlying physical 
principle of the nonlinear evolution is the tendency for 
the growing fluctuations to compensate on the average 
the pressure anisotropics generated by the turbulence. 

This mechanism of making small-scale magnetic fields 
is distinct from the fluctuation dynamo, which exponenti- 
ates the magnetic energy at the turbulent stretching rate 
(~ 70, much slower than the plasma instabilities) and 
produces long filamentary folded structures, so the paral- 
lel correlation length of the field remains macroscopically 
large (~ outer scale) [13j, in contrast to the instability- 
produced wrinkles with parallel scales possibly as small 
as the ion gyroscale. How the dynamo operates in the 
presence of the instabilities 3 is a subject of an ongoing 
investigation motivated by the fundamental problem of 
the origin of cosmic magnetism in general and of mag- 
netic fields in galaxy clusters in particular. 

To illustrate the potentially dramatic effect of fire- 
hose and mirror fluctuations on the transport proper- 
ties of magnetized turbulent plasmas, consider the ICM 
thermal conduction problem. The standard estimates of 
the electron thermal conductivity in a tangled magnetic 
wtheAmfp if A m f p < l B (collisional), 
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field are 

K e ~ wth e A m fp^/LRR if l B < A mfp < L RR (semicol- 
lisional), and n e ~ «the^B if A m f p 3> Lrr (collision- 
less), where Ib is the (parallel) correlation length of the 
magnetic field and Lrr = IbIuQb/ Pe) is the Rechester- 
Rosenbluth length. In most MHD models [2l[ (including 
the fluctuation dynamo [Hj]), Ib is macroscopic and all 
three estimates yield an effectively isothermal ICM (ex- 
cept at macroscopic scales). However, if magnetic wrin- 
kles with 5B/B ~ 1 develop at scales ~ pi, we have 
Ib ~ Pi and Lrr ~ pih\(pi/p e ) < A mfp , so n e ~ v t hepi- 
For our fiducial Hydra A parameters, this is 10 10 times 
smaller than the collisional value, so there is effectively 
no thermal conduction on macroscopic scales. The ICM 
viscosity is similarly reduced, from t>thiA m f p to VthiPi be- 
cause with I b ~ Pi, the effective ion mean free path is 
~ pi . Curiously, in stronger-field regions where 2 //3, > A 
and the instabilities are suppressed, the transport is more 
effective: the thermal conductivity and viscosity remain 
large (although highly anisotropic). 

Due to spatial resolution constraints, the firehose and 
mirror structures are not directly detectable in clusters, 
but the huge changes in the transport coefficients that 
they may cause will have a potentially predictable effect 
on observable large-scale fields and flows [H, [f| . More di- 
rect information is available from satellite measurements 
in space plasmas. Mirror structures with 5B/B ~ 1 have, 
indeed, been found [22| and there is strong evidence that 
the directly measured temperature anisotropics match 
the firehose and mirror marginal stability conditions [23| . 
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